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ABSTRACT 

The enigma source, RX J1856. 5-3754, is one of the so-called dim thermal 
neutron stars. Two puzzles of RXJ1856. 5-3754 exist: (1) the observational X-ray 
spectrum is completely featureless; (2) the UV-optical intensity is about seven 
times larger than that given by the continuation of the blackbody model yielded 
by the X-ray data. Both the puzzles would not exist anymore if RX J1856. 5-3754 
is a low mass bare strange quark star, which is in a propeller phase with a low 
accretion rate. A boundary layer of RX J1856. 5-3754 is suggested and modelled, 
from which the UV-optical emission is radiated. Free-free absorption dominates 
the opacity of the boundary layer, which results in the opacity to be high in 
UV-optical but low in X-ray bands. The star's magnetic field, spin period, as 
well as the accretion rate are constrained by observations. 

Subject headings: dense matter — elementary particles — pulsars: general - 
stars: neutron 



1. Introduction 

Quarks (and leptons) are fundamental fermions in the standard model of particle physics, 
the underlying theory of the interaction between which is believed to be quantum chromo- 
dynamics (QCD). Quark matter (or quark-gluon plasma) is expected in QCD, but is not 
uncovered with certainty. While physicists are trying hard to find quark matter by ground- 
based colliders, astronomers are researching on the existence and/or consequence of quark 
matter in the sky. Quark (matter) stars may form by supernova explosions, which are sug- 
gested to manifest as pulsar-like compact stars. RX J1856. 5-3754 (hereafter RX J1856) is one 
of the .ROS^T-discovered so-called dim thermal neutron stars. It was thought to be actually 
a quark star (Drake, J. et al. 2002; Xu 2002). But this idea was soon questioned by Walter 
& Lattimer (2002) if a two-blackbody model is assumed in order to fit the observations by 
Chandra, EUVE, and HST. 

Among the isolated compact objects, RX J1856 is the mostly studied one because it is 
the brightest. The parallax distance of this source is only 117 ± 12 pc (Walter & Lattimer 
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2002) . Its optical counterpart, with V~26 mag, was found (Walter & Matthews 1997), but no 
radio counterpart has been observed. A 505-ks Chandra observation revealed a featureless 
spectrum which can be well fitted by a blackbody, with apparent radius Roo= 4.4 km 
and temperature T 00 = 63 eV (Burwitz et al. 2003), and no X-ray modulation being found. 
Ransom et al. (2002) setted an upper limit on the pulsed fraction of ~4.5% (99% confidence) 
for frequency < 50 Hz and frequency derivative —5 x 10 10 Hz s _1 < / < Hz s _1 , whereas 
Burwitz et al. (2003) obtained an upper limit of 1.3% (2a confidence) on the pulsed fraction 
in the frequency range (10~ 3 — 50) Hz using XMM-Newton data. 

Various efforts have been made to explain the observations in the regime of normal 
neutron star, such as a two-component blackbody (Triimper et al. 2003), a neutron star with 
reflective surface (Burwitz et al. 2003), a magnetar with high kick- velocity (Mori & Ruderman 

2003) , a naked neutron star (Turolla et al. 2004), and a surface with strong magnetic field 
(van Adelsberg et al. 2005; Perez- Azorin, Miralles & Pons 2005). But these models are far 
away from fitting the real data with reasonable emissivity. Summarily, several difficulties 
are in the normal neutron star models. (1) Featureless spectrum could be hard to reproduce 
by a neutron star with an atmosphere of normal matter. Though strong magnetic field may 
smear out some spectral features (Lai 2001), the consequent large spin down luminosity is 
not observed. (2) Special geometry is needed to explain the no-modulation observation - 
the pulsar should be aligned or we are situated at the star's polar direction. (3) A normal 
neutron star of radius 17 km in the two-component model requires a very low-mass about 
O.4M . The mechanism to form such a normal neutron star is still unknown. 

Alternatively, the X-ray observation alone may be understood by assuming RX J1856 to 
be a low-mass bare strange quark star, and Zhang, Xu & Zhang (2004) fitted well the X-ray 
data with a phenomenological spectral emissivity in a solid quark star model. However, the 
UV-optical observations revealed a seven times brighter source comparing to that derived 
from X-ray observation. What's missed here? To overcome these difficulties, we present a 
model under low-mass quark star regime. A boundary layer around RX J 1856 is proposed 
in this paper, which is optically thick for UV-optical radiation but is optically thin for X- 
ray radiation. We can obtain consistency with the observational data, as well as strong 
constrains on the star's magnetic field strength and spin period, through this approach. 



2. The model 

We consider the star to be a bare quark star which would be indicated by the featureless 
spectrum from the 505-ks Chandra data (Xu 2002). Since a bare quark star have quark 
surface instead of atmosphere, it could reproduce a Planck-like spectrum, rather than a 
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spectrum with atomic lines (Xu 2003). Additionally, the small apparent radius (i?^ = 4.4 
km) observed in X-ray band could be another hint. For these reasons, a bare quark star 
hypothesis might be reasonable. The following discussions are under this hypothesis. 

The object, RX J1856, in our model consists of two components: (1) a central bare quark 
star which radiates X-ray photons; (2) a boundary layer at the magnetic radius r m (Alfven 
radius) with a quasi thermal spectrum. There are also two free parameters in our model: (1) 
/i m , the magnetic moment per unit mass of the quark star; (2) M, the accretion rate. Strong 
constrains on these two parameters will be given by observations in this model. Accreted 
matter should be stopped at the magnetosphere and form a boundary layer. According to 
following calculations, the boundary layer is optically thin for the central X-ray emission but 
thick for UV-optical emission because the cross section is much smaller for X-ray photons 
than UV-optical ones. By this model, we can fit the observations in UV-optical and X-ray 
bands. The details are below. 



2.1. The star 

For a low-mass quark star, the internal density is almost homogenous (Alcock et al. 
1986). The star's mass and density can be well approximated as 

4 

M ~ -iiR 3 p and p ~ 4B, (1) 
o 

respectively, where the bag constant is suggested to be B = (60 ~ 100) MeV fm~ 3 = 
1.07 ~ 1.96 x 10 14 g cm" 3 . A median bag constant, B = 1.5 x 10 14 g cm" 3 , R^ 4.4 
km and T OQ = 63 eV are adopted in this paper. Considering general relativistic effects 
(Roo = R/(l - RJR) 1 / 2 ), we have R = 4.3 km, M = 0.097M©, and R g = 2GM/c 2 = 0.29 
km. Since M is much less than 1M , the general relativistic effects are rather small. The 
magnetic moment \i of the star could be expressed as (Xu 2005) 

H = ^BR 3 = M/v. (2) 

This is a relatively free parameter because it spans a large range, (10~ 4 ~ 10~ 6 ) G cm 3 g _1 
(Xu 2005). We use /i m as a free parameter which will be constrained strongly. 

2.2. The boundary layer 

Since we have never observed any accretion powers, such as strong X-ray emission and 
X-ray burst, the accretion rate of the star should be very low. Hence we consider ADAF 
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(advection dominated accretion flow) model for the accretion (e.g., Narayan et al. 1998), 
which means M <~ 0.01M E dd, where M Edd is Eddington accretion rate of RX J1856, 

M Edd = 4vrm pCj R/(T T ~ 4.2 x 10 17 g s" 1 , (3) 

with m p the proton mass and <tt the Thomson cross section. The star should be in a 
propeller phase, and most of the accretion matter can not be accreted onto the the star's 
surface (otherwise the accretion induced X-ray luminosity will be much higher than that we 
observed) . 

By the accretion flow, a quasi spherical layer may form at the magnetospheric radius, 

R2 d6 

r ~ ( ) 2/7 (A) 

which is derived by equaling kinematic energy density of free-fall matter to magnetic energy 
density. At the magnetospheric radius, r m , accreted matter will decelerate and pile up to 
form a boundary layer. The density at the boundary layer should be high enough to make 
it an optically thick region. The emission spectrum could then be blackbody-like. 



2.3. Constrain from temperature 

The UV-optical observation of RX J1856 is in the Rayleigh- Jeans tail of the radiation 
from the boundary layer. The brightness of a back body for hv <C kT is B v ~ 2ku 2 T rn Jc 2 , 
where k is Boltzman constant. The observed flux per unit frequency should then be 

F„ = nByjd 2 ~ 2irk^T rin r 2 m , (5) 

where d is the the star's distance to the earth. Observationally, one has 

T rin r 2 n = const. = 9.5 x 10 3 eV km 2 . (6) 

In order to be consistent with observation, the boundary layer's temperature should be 
neither too high (not to affect the X-ray spectrum) nor too low (to keep the Rayleigh- Jeans 
slope), i.e. (Burwitz et al. 2003), 

4 eV < T rm < 33 eV. (7) 

From Eqs.(2), (4), (6), and (7), one comes to 

AM 2 n 2 

17 km < ( ^ N2/7 < 49 k (g) 
MV2GM 

where only /i m and M are free parameters if the star's radius (and thus the mass through 
Eq.(l)) is determined observationally. These two inequations of Eq.(8) give two dashed lines 
in M-/i m diagram and constrain the allowed region to a small belt (see Fig. 1). 
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2.4. Constrain from optical depth 

To reproduce the observed spectrum in UV-optical bands, the boundary layer should 
be optically thick. We assume that the accreted matter's radial velocity decreases from free 
fall velocity vg to zero at the boundary layer. From mass continuity pv = M/47rr^, the 
density p should increase inward. The magnetic pressure at radius r is V(r) = B(r) 2 / (8ir) = 
p? j (27rr 6 ). Let's Consider a cubic mass with border a and density p near the boundary layer 
of the mass flow. The cube would feel a force by the magnetic field 

F~g«". (9) 

The pressure gradient dP jdr oc r ~ 7 is proximately a constant if the the layer is thin. With 
the conservation of radial momentum, — Fdt = d(mv), and the definition of differential 
displacement, dx = —fdt, we have Fdx ~ mvdv ~ a 3 pvdv. One could then estimate the 
x-value of the layer (x ~ at the bottom where v ~ and p = p max ) to be, 

M dV 

x~v- — T l— — oc v. (10) 

^ixr 2 J Or v ' 

We obtain thus p oc 1/x since pv is a constant. It would not be easy to estimate p max since 
the coupling between the accretion disk and the star's magnetosphere is difficult to model. 
We have thus another free parameter, p max , in our calculations. Note that v is the radial 
velocity. The total velocity does not decrease so much because the velocity changes direction 
and most of the accreted matter will be propelled out. 

At temperature being higher than 4 eV, one can calculate from Saha equation that nearly 
all hydrogen atoms is ionized. The opacity is produced mainly by Thomson scattering, bond- 
free and free-free absorptions. In the following, we show that free-free absorption dominates. 

The Thomson cross section <tt is a constant, and the optical depth is very small if 
only Thomson scattering is considered. Free-free absorption is sensitive to density, so it is 
dominated in the high density region and may result in an optically thick boundary layer. 
Since free-free absorption coefficient (Rybicki & Lightman 1979, p. 162) 

kb = 3.7 x l0 8 T-^ 2 Z 2 n e niU- 3 (l - exp(~))g s , (11) 

where n c and n\ are electron and ion number density, Z ~ 1 is ion charge and gg ~ 1 is 
Gaunt factor, optical depth differs by a factor of ~ 10 6 between UV-optical and X-ray bands. 
Therefore it is possible that the boundary layer would be always optically thin for the central 
X-ray emission. Detail calculations confirm this suggestion (see Fig. 2). 
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Bond-free and bond-bond absorptions should also be considered because these processes 
may induce features in the spectrum which are not observed yet. An estimation goes as 
following. For oxygen, the most abundant element after hydrogen and helium, the typical 
bond-free cross section ~ 10~ 19 cm 2 in 0.1-1 keV band (Reilman & Manson 1979). Taking 
into account its abundance ~ 10~ 3 (number relative to hydrogen), the effective cross section 
is a"bb-bf ~ 1CT 22 cm 2 . The optical depth induced by these absorptions could be smaller 
than 1 though Obb-bf is about ~ 10 2 times that of electrons. The model proposed may then 
still work. 



3. The results 

3.1. The radiation efficiency in the layer 

Since the star is in a propeller phase, most of the accreted matter will be expelled out 
finally. In this way only a small part of the accretion energy may be converted to radiation. 
Here we use / ra d to represent the transformation rate. We have then 

GMM 4 2 

ikd = vT n BL r m , (12) 

where Qbl is the emission angle of the boundary layer which is An if the boundary layer is 
completely spherical. Since / ra( j and ^bl are degenerate, we adopt simply Qbl = 47r in the 
following. 

Combining with Eqs.(2), (4), (6), and (12), one may obtain the value of a new free 
parameter / ra d, the dependence of which on other two free parameters is 

/ rad ~ 9.24 x 1(T 32 ^— . (13) 

/im 

Different / rad value gives different lines on M-pL m diagram (see Fig. 1). The allowed range 
of /rad should be about 10~ 2 ~ 10~ 4 . 



3.2. The star's spin period and polar magnetic field 

A propeller phase requires r co < r m , where r co = (GMP 2 /(47r 2 )) 1 / 3 ~ 6.88 x 10 7 P 2 / 3 cm 
is the corotating radius. We can thus constrain the star's spin period to be < 20 ms from 
this requirement. On the P — P diagram, a death line may simply be represented by an 
equal voltage line (Zhang 2003). We therefore constrain P and P of RX J1856 via setting 
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Fig. 1. — The magnetic moment per unit mass /i m (G cm 3 g _1 ) vs. the accretion rate 
dM/dt (same to M, in units of Eddington accretion rate). The dashed line are constrains 
from Eq.(8). Lines labelled by different r(10eV), the optical depth of the layer for his = 10 
eV, are drawn via Eq.(ll). The dotted line is for r(O.lkeV) = 0.1 with hv = 0.1 keV, through 
Eq.(ll) too. Lines denoted by / ra( j are for different radiation efficiencies, from Eq.(13). The 
allowed parameter space is limited by four lines: the two dashed lines, the dotted line, and 
the line denoted by r(10eV) = 1. It is evident that 8 x 10~ 9 < /i m (G cm 3 g" 1 ) < 8 x 10~ 8 , 
10~ 4 ~ /rad ^ 10 -2 , and an accretion rate 10" 3 < M/M Edd < 10~ 2 which is consistent with 
the ADAF assumption. 
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Fig. 2. — Plot of optical depth r vs. photon energy E. The bottom density, p max , is 
considered as a free parameter. Two hatched regions are for UV-optical and X-ray ob- 
servations, respectively. The dotted line is for r = 1, beyond which the layer could be 
optically thick. Thick lines are for free-free absorption, while thin lines for Thomson scat- 
tering. Solid, dash-dotted and dashed lines are for different p max - It is found that the value 
Pmax = (10 5 ~ 10 6 )p(r m ) in order to reproduce a boundary layer to be optically thick in UV- 
optical but thin in X-ray bands, where p(r m ) = M / (47rr m Vff(r m )). We choose M = lCT 3 MEdd 
and p m = 2 x 10~ 8 G cm 3 g _1 for indications here. 
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the death voltage to be 10 12 V in order to produce approximately the dash-dotted death line 
for normal pulsars in Fig.3. The death line would be described by 

$ ~ 6.6 x 10 12 B 12 P- 2 RI V ~ 10 12 V, (14) 

where B u is the magnetic field strength in units of 10 12 G and R 6 is star radius in units of 
10 6 cm. The line for RX J1856 is then 0.5B 12 P~ 2 ~ 1. Note that since low-mass quark stars 
could have small radii, their death line would be different from that of the normal pulsars. 
A constrain for magnetic polar field (Xu 2005), B = 32irBfi m /3, is obtained so as to make 
the star radio quiet, 4 x 10 7 G <B<4x 10 8 G, from the limit 8 x 10" 9 < /i m (G cm 3 g _1 ) < 
8 x 10~ 8 . We could thus obtain 10~ 22 < P < 10 -20 in Fig. 3 if a spindown torque by pure 
magnetodipole radiation is assumed. 



4. Conclusions and discussions 

We propose that RX J1856. 5-3754 is a low-mass quark star, which is in a propeller 
phase. A boundary layer of RX J1856 is suggested and modelled in order to understand 
the UV-optical excess observed. In order that the layer could be optically thin in X-ray but 
thick in UV-optical bands, some constrains for RX J 1856 are obtained: the polar magnetic 
field 4 x 10 7 G <B<4x 10 8 G, the spin period P < 20 ms, and an accretion rate 10 -3 < 
M/M E dd < 10 -2 , if the the star's radius and mass are 4.3 km and ~ O.1M , respectively. 
Additionally, we find that the radiation efficiency of accretion matter in the layer could be 
between 10~ 2 and 10~ 4 , and that the density at the layer's bottom is about 10 5 ~ 6 times 
that at the top. No magnetospheric activity exists on the star since it would be under its 
death line, which could be the reason that neither radio emission nor X-ray modulation is 
observed from RX J1856. The star could have a weak magnetic field to be similar to that of 
the millisecond pulsars. 

In order to reproduce the UV-optical spectrum, an optically thick layer is considered 
only, which results naturally a slope of the Rayleigh- Jeans tail. We are still not sure whether 
the UV-optical spectrum observed could also be understood by an optically thin boundary 
layer. The difficulty in the later case could be to calculate the emissivity of matter with strong 
magnetic field. Even in the first case of an optically thick layer, the radiation efficiency is 
still hard to calculate due to few knowledge of the interaction between charged particles and 
magnetic field. We obtain a low radiation efficiency, 10 -4 2 , in our constrains, which could 
be reasonable since most of the kinematic energy is take out during a propeller process. The 
emission from the outer part of the ADAF disk is not taken into account since the luminosity 
should be very low, with emission in low-energy band. 
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Fig. 3. — The P — P diagram and the possible position of RX J1856. RX J1856 should 
be in the left of the solid vertical line because of P < 20 ms. It should also be in-between 
two dotted lines through 4 x 10 7 < BG < 4 x 10 8 G. The dash-dotted and dashed lines 
are approximately death lines for normal pulsar with radius R = 10 km and for RX J 1856 
with R = 4.3 km, respectively. It is suggested that RX J 1856 could be in the triangle 
(i.e., the "allowed region") if the star spins down due to pure magnetodipole radiation. No 
magnetospheric activity exists on the star because of a low potential drop in the open field 
line region. The pulsar data are from ATNF Pulsar Catalogue (Manchester et al. 2005, 
http : // www. atnf . csiro . au / research / pulsar / psrcat ) . 
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The bag constant (or effectively, the average density of the star AB) is chosen to be 
B — 1.5 x 10 14 g cm" 3 in our calculations. However, the allowed region in P — P diagram 
will be larger if B is released from 1.07xl0 14 g cm" 3 to 1.96xl0 14 g cm -3 . The model 
proposed in this paper would thus be a flexible one. 

How to check the model proposed? A direct way might be to find the star's spin period 
in radio and/or X-ray data. This model would be rule out if the spin period is much larger 
than 20ms. Photons in low-energy bands (e.g., infrared or submillimeter emission) would also 
be radiated from the ADAF disk and the boundary layer. The luminosity of this emission 
could be high enough if more powerful facilities are offered in the future. 
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